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ABSTRACT 
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A procedure  for  the  numerical  approx imatian  of  the  Cauchy  problem  for 
the  following  linear  parabolic  partial  differential  equation  is  defined: 

ufc  - (p(x)ux)x  + q (x)  u = 0,  0 < x < 1 , 0 < t <_  T;  u (0,t)  = f^t)  , 

0 < t <_  T;  U (l,t)  = f2(t),  0 < t <_  T;  p(0)u  (0,t)  = g(t), 

0 < t <_  t <_  T;  | u (X,t)  | <_  M , 0 < x < 1,  0 <_  t <_  T . 

The  procedure  involves  Galerkin  type  numerical  methods  for  related  parabolic 
initial  boundary -value  problems  and  a linear  programming  problem.  Explicit 
a priori  error  estimates  are  presented  for  the  entire  discrete  procedure 
when  the  data  f^,f2,  9 are  known  only  approximately. 
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SIGNIFICANCE  AND  EXPLANATION 

In  many  physical  problems  in  heat  conduction,  it  is  impossible  to  obtain 
an  initial  temperature  distribution  within  a material.  In  many  of  these 
cases,  in  order  to  obtain  approximations  of  the  temperature  within  the  body, 
one  must  rely  entirely  upon  data  which  can  be  measured  at  the  boundary. 

An  additional  problem  is  that  these  boundary  data  are  only  accurate  to  within 
some  prescribed  measurement  errors. 

The  purpose  of  this  paper  is  to  define  a procedure  for  numerically 
approximating  the  solution  of  one  such  heat  flow  problem  and  to  present 
explicit  error  estimates  for  the  numerical  procedure.  A priori  error 
estimates  are  presented  when  the  data  are  known  only  approximately. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
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NUMERICAL  APPROXIMATION  OF  A CAUCHY  PROBLEM 
FOR  A PARABOLIC  PARTIAL  DIFFERENTIAL  EQUATION 


Richard  E.  Ewing  and  Richard  S.  Falk 


1.  Introduction. 


Consider  the  numerical  approximation  of  the  solution  of  the  following  Cauchy  problem 
for  a linear  parabolic  partial  differential  equation. 

Problem  (P) : Find  a function  u = u(x,t)  satisfying 

u = 0 , 0 < x < 1 , 0 < t £ T , 

0 < t <_  T , 

0 < t < T, 


a) 

3u  3 , , , 3u  . 

at“  37  (p(x)  37  > + 

b) 

u(0,t)  = f (t) 

(1.1) 

c) 

u(l,t)  = f2  (t) , 

d) 

p(°)  |^-  W.t)  = g(t) , 

e) 

|u(x,t)  I <_  M , 

where 

the 

data  f^,  and  g 

a) 

lfl-f?[O.T]  ^£0  ' 

(1.2) 

b) 

11  f2  _ f2!l  (0 ,T]  - £0  ' 

c) 

" 9 - 9*"tt0,T]  1£0  , 

with 

Eo  5 

■ 0 and  where  for  any 

(1.3) 

11  f"  [a,b]  = ^ 

a < t < b 

0 < X < 1 , 


0 ' to—  1 1 T » 
0 1 1 T * 
as  fj,f*,  and 


|f(t)|- 

We  assume  that  the  following  hypotheses  are  satisfied: 

(HI)  and  q are  such  that  a classical  solution  u to  (1.1)  exists. 

(H2)  The  functions  p,  p'  and  q are  uniformly  Holder  continuous  in  0 _<  x <_  1 
and  satisfy 

a)  0 < p*  < p(x)  <_  p*  , 

b)  0<q,  < q(x)  < q*  , 

c)  |p'  (x)  [ p ' * , 

(H3)  f^  and  f are  continuously  differentiable  and  a constant  exists  such 
that 

11  f 1H  [0,T]  + 11  fX11  [0,T]  + "V  [0,T]  + " f2M  [0,T1  + llg1l[t0>T]  ‘ 

Sponsored  by  the  United  States  Army  under  Con tract Numbers  DAAG29-75-C-0024  and  DAAG29-78 - 
G-0161.  This  material  is  based  upon  work,  supported  by  the  National  Science  Foundation 
under  Grant  Numbers  MCS78-09525  and  MCS78-02  737. 
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(H4)  For  the  Sturm-Liouville  problem  with  eigenvalues  and  corresponding  normalized 
eigenfunctions,  satisfying 

a)  (p^')'-qy’  + X = 0,  0 < x < 1 , 

n n n n 

b)  yp  (0)  = p (1)  = 0 , 

n n 

there  exists  a d > 0 such  that 


d = inf  (X 


n+1 


X ). 
n 


(H5)  The  approximations  (1.2)  hold.  In  addition 

a)  11  fi  + T'ho.X]  - 


b) 


If'  - f * • II  < c . 

2 2 ' [0,T]  - 0 


(H6)  If  M > 0 is  the  constant  given  in  (l.l.e),  we  have 
3)  " fl"  [0,T]  — M ' 


b) 


1 f2l!  [0,T]  - M 


C)  llfi"[0,Tl  ' 

d)  "^"[O.T]  • 

Without  the  extra  assumption  (l.l.e),  the  Cauchy  problem  (1.1. a-d)  is  not  well- 
posed  in  the  sense  of  Hadamard  [3,4,5,6,7,8,13]  since  the  solution  does  not  depend  contin- 
uously upon  the  data.  However,  for  the  problem  as  stated  with  hyptheses  (H1)-(H6) 
satisfied,  a continuous  dependence  result  was  obtained  in  (6]  . 

Many  studies  of  the  continuous  dependence  of  the  Cauchy  problem  for  various  para- 
bolic problems  have  appeared  in  the  literature  [3-8.  13].  Other  formulations  in  terms  of 
control  problems  have  also  appeared  in  the  control  theory  literature  [11,  12,  14).  In  [7] 
Ginsberg  considered  numerical  approximation  of  the  Cauchy  problem  for  the  heat  equation 
ut  - u^x  with  g(t)  = 0 in  (l.l.e)  by  expanding  the  data  in  Fourier  series  and  estimating 
Fourier  coefficients.  In  [3,4]  Cannon  and  Douglas  outlined  numerical  procedures  for  various 
Cauchy  problems  for  the  heat  equation  by  reducing  the  probl  ms  to  mathematical  programming 
techniques.  In  [5]  Cannon  and  one  of  the  authors  presented  a direct  numerical  method  for  a 
slightly  different  Cauchy  problem  for  the  heat  equation  in  which  a Taylor  series  expansion 
for  the  data  is  numerically  approximated.  In  [6]  a numerical  scheme  requiring  numerical 
approximation  of  several  unknown  eigenvalues  and  eigenfunctions  was  presented  without 
explicit  error  estimates  for  the  approximations.  In  this  paper  the  numerical  schemes 
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involve  only  solution  of  linear  parabolic  initial  boundary-value  problems  and  a simple 
linear  programming  problem.  More  importantly,  explicit  a priori  error  estimates  for  the 
entire  procedure  are  presented. 
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In  Section  2,  basic  notation  is  presented  ana  Problem  (P)  is  reduced  by  linearity 
into  two  simple  initial  boundary-value  problems  and  an  optimization  problem.  In  Section  3, 
Galerkin  type  numerical  schemes  are  defined  for  the  initial  boundary -value  problems  and 
a linear  programming  problem  is  formulated  to  solve  the  optimization  problem.  Several 
basic  lemmas  needed  to  prove  the  main  result  are  stated  in  Section  4.  Then  a priori  error 
estimates  for  obtaining  approximations  to  (1.1)  with  approximate  data  satisfying  (1.2)  are 


2 . Preliminaries. 


We  shall  first  define  some  of  the  notations  used  for  various  norms  throughout  the 
paper.  Recall  that  in  Section  1 we  used  the  notation  that,  for  any  function  f = f(t)  , 


II  fll 


la.b] 


sup  | f (t)  | 
a < t < b 


For  functions  ij/  = i|j(x)  defined  on  (0,1),  we  shall  denote  by  |i|i|  the  norm  II 


L (0,1) 


and  by  W 


,m,“ 


(m  a positive  integer)  the  usual  Sobolev  space  of  functions  with  norm 
m 


l 

j=0 


aJtp 


3x 


s , 2 


For  real  s , we  further  denote  by  H the  Sobolev  space  W ' of  real-valued  functions 

2 

defined  on  (0,1)  and  by  II  i|ill  its  corresponding  norm.  We  note  that  L (0,1)  will  be 


denoted  by  H°  and  II 1 1/ 


II  by  II  i|)  II  . For  definitions  of  the  other  spaces,  we  refer 

L2(0,1) 


the  reader  to  [101  . 

Also,  for  X a normed  space  with  norm 


and  u:[a,b]  -*•  X , we  define 


Hull 


L (a,b;X)  a 


s Jbiiu(-,t)ii; 


dt 


and 


Hull  ^ = sup  llu(*  ,t)ll  . 

L (a,b;X)  a <_  t <_  b 

Finally,  for  convenience,  we  define  a bilinear  form 

a (u,v)  5 (p(-)ux,vx)  + (q(*)u,v)  , 

2 

where  (•,•)  denotes  the  L (0,1)  inner  product. 

We  shall  next  present  a reformulation  of  Problem  (P)  on  which  our  approximation 
scheme  will  be  based.  We  first  define  a function  w satisfying  the  initial  boundary 
value  problem 

a)  If  - h (p  S’ + qw=  F(x'fc)'  0 * x < 1(  0<t^T' 

b)  w (0,t ) = 0 , 0 < t <_  T , 

(2.1) 

c)  w ( 1 ,t ) = 0 , 0 < t <_  T , 

d)  w(x,0)  = 0 , 0 < x < 1 , 

where 


F (x,t ) = 


(2.2) 


-(l-x)q(t)  - xf  2 < t ) + p’(x)[f2(t)  - fx  (t)  1 - q[(l-x)f1(t)  + xf  2 (t)  1 . 


2 ii> 

Also,  for  each  ^ e L (0,1),  we  define  a function  z satisfying  the  initial  boundary 


value  problem: 

a) 

b)  zv(0,t)  = 0 , 

(2.3) 


Sr*  3 , Sr*  <|i 

if  - *r  (p  it  > + * 


0 , 


0 < x < 1 , 

0 < t <_  T , 

0 < t <_  T , 

0 < X < 1 . 


0 < t < T , 


c)  zMl,t>  = 0 , 

d)  z^(x,0)  = i|/  , 

Using  M from  (l.l.d),  let 

K 2 {*  € C°(0,H:  1*1^  <_M  and  V>  (0)  = f.^0),  <P  (1)  = ^(0)}. 

We  can  now  reformulate  Problem  (P>  using  (2.1)-(2.3)  and  linearity  of  the  operator  in  (1.1) 
as  follows: 

Find  <f  (corresponding  to  u(x,0)  from  (1.1))  with  £ K such  that 


P (0) 


3z 


V-i 


3w 


(2.4) 


3x  (O.t)  = g (t)  - p (0)  (0,t)  - p(0)[f2(t)  - fx(t)] 

5 G(t) 


r (x)  5 (l-x)f1(0)  + xf2(0)  . 


where 
(2.5) 

To  see  that  this  reformulation  is  equivalent  to  Problem  (P),  observe  that,  using 
linearity,  the  function 

= W + Z^  + (l-x)f^(t)  + Xf2(t) 


satisfies 

a) 


3u*  3 , 3u*  + * 

if  ■ 3x  (p  ~nr)  ^ ■ 0 

j, 


0 < x < 1, 

0 < t <_  T, 

0 < t < T , 


0 < t < T, 


0 < X < 1 , 


b)  uv(0,t ) = f^ft)  , 

(2.6)  c)  u*(l,t)  •=  f2  (t)  , 

d)  u^fx.O)  = <l>  + (l-xlf^O)  + xf 2 (0)  , 

e)  p(0)u*(0,t)  = p(0)wx(0,t)  + p(0)z*(0,t>  + p (0)  t f 2 (t)  - f ^ (t ) ] , tQ  < t < T. 

J)  — jr 

Choosing  = y>-r,  we  see  that  u solves  problem  (P)  provided  we  can  find  a v>  satis- 
fying (2.4).  (We  note  that  since  | f x (t > | <_  M and  | f 2 (t ) | <_  M,  0 <_  t <_T,  the  maximum 
principle  implies  that  the  condition  }^>  | ^ M is  equivalent  to  requiring  that 

( u (x * t ) | < M,  for  0 < x < 1 , 0 < t < T.) 
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Using  this  notation,  we  can  now  state  a form  of  the  continuous  dependence  result 
(proved  in  [6]  ) which  we  shall  need  later  in  the  derivation  of  the  error  estimates. 
Lenina  1 : Under  hypotheses  (H1)-(H6),  there  exist  computable  constants  and  y 

(0  < y < 1)  such  that  for  all  ip  with  < M , 

(2-7)  ll2*"  - 1,~  iCl  llzx(0'-)llIt  T1  ' 

l b 

where  z is  the  solution  of  (2.3). 
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3.  Description  of  tla-  N ume  r u:a  l App  r ox  ima  t ion  s . 

Tn  this  section  we  consider  the  problem  of  numerically  approximating  the  solution  of 

(1.1)  (or  equivalently,  its  reformulation  described  in  Section  2)  subject  to  the  restriction 

(1.2) .  The  restrict icn  (1.2)  comes  from  the  fact  that  data  measurement  error  is,  in  general, 

accurate  only  to  within  some  measurement  tolerance  . We  denote  by  F*  the  function 

obtained  by  replacing  f ^ and  f ^ , in  (2.2),  by  f*  and  f*  respectively.  We  further 
define  w*  as  the  solution  of  (2.1)  with  F replaced  by  the  F*  defined  above. 


(3.1)  r * (x)  - (1-x) f * (0)  + xf* (0) 


(3.2)  G*  (t ) g*  (t ) - p (0)  | (0,t)  - p(0)[f*(t)  - f * (t)  ] . 

d X Z L 

We  note  here  that  from  (H5),(H6),  (2.5),  and  (3.1)  we  have 

a)  | r | < M 

(3.3)  b)  I r*  i < M 

c)  II  r-r * II  : ! r-r*  1 

0 — ' 'oo 

<_  max  { (1-x)  | f x (0)  - f*(0)  | + x|f2(0)  - f * (0)  | } 

0 < x < 1 

- eo  ' 

Since  the  data  is  only  known  approximately  as  described  above,  we  now  define  an 
approximation  sc  lie  me  based  on  a finite  dimensional  analogue  of  equation  (2.4)  with  G 
replaced  by  G*  from  (3.2).  We  first  describe  schemes  for  obtaining  approximations  to 

i if 

w and  zr  (assuming  ip  is  known). 

o 

Let  sj  denote  the  space  of  continuous  piecewise  polynomials  of  degree  6-1 
defined  on  a uniform  mesh  of  width  h on  10,1)  and  vanishing  at  x = 0 and  x = 1. 

We  shall  consider  a family  of  such  spaces  for  0 < h < 1 . We  assume  each  space  in  this 

g 

family  satisfies  the  following  so-called  "inverse  assumptions":  u e implies  that  for 


some  constant  , 

a)  |u|  „ < C2h_1|uL. 

(3.4)  L'  - 

b)  1 u | < C h II  tall  _ - 

i i oo  — 2 0 

a 

We  also  assume  that  our  test  spaces  satisfy  the  following  approximation  assumptions: 

8 8 
if  u c S,  and  I,u  is  the  interpolate  of  u in  S,  , then  for  some  constant  C , 
h h n 3 
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a) 

(3.5) 

b) 

llu  - 

iv 

Let 

k > 0 

We  shall  present  a 

Define  W : 

C-l 

o 

II 

ft 

Z.ull  < C,  h Hull 
n 0 — 3 


s > 0 , 


< C,  u 


1»“ 


:<t"). 


We  shall  present  a Crank-Nicolson-Galerkin  approximation  for  w , the  solution  of  (2.1)- (2.2), 

_4 


0 1' 


N 


by 


(3.6) 


w — w w + w 

; n-+1-k  ■ n , v)  + a ( n+17  n , vl=  (F  (.,(»+  j)k),V) 


for  all  V € and  n = 0,1 , . . . ,Nt~1  , with  WQ  = 0 , where  a(«,«)  is  defined  in 

Section  2.  We  similarly  define  VP  to  be  the  analogous  approximation  to  w*(nk)  given 
by  (2.1)- (2.2)  with  F replaced  by  F*. 

The  scheme  defined  in  (3.6)  is  known  to  have  a time-truncation  error  of  the  order 
2 2 

k . We  shall  use  another  0(k  ) time-stepping  method  with  better  stability  properties  but 

greater  work  estimates  to  approximate  z , the  solution  of  (2.3)  with  ip  assumed  known. 

Define  Z*:  {0  = t ,t  , ,t  = T)  -v  by 

u 1 n 

T 


(3.7) 


n+a  n „ 'i  . _ ,„\p 


ak 


~ 'V  ) + 3 (Zn+a'v)  = ° ' 


z r — z T 

f n+1  n i . ,„tj i 

1 S 'vJ  + a(Zn+l'V) 


a(z^+  ,v)  , 
a n+a 


with 


(3.8)  (z0'v)  = (,(''v)  ' 

2 /2 

for  all  v £ S,  where  a = 1 - — . 

k 2 

2 

We  note  that  since  each  of  the  time  stepping  schemes  defined  above  have  0(k  ) 
time-discretization  error  but  different  spacial  orders  of  approximation,  we  shall  use  the 
time  step  k to  tie  the  two  approximations  together.  Thus  k will  be  the  same  in  each 
of  (3.6)  and  (3.7).  We  shall  then  see  that  in  order  to  balance  the  temporal  and  spatial 

discretization  errors  in  each  problem  separately,  we  shall  let  h = k in  the  definition 
1/2 

of  Z and  h^  = k in  the  definition  of  W . 

Let  Ng  = ItQ/kl+  1 where  IrJ,  for  re®,  is  the  greatest  integer  less  than  r . 
Using  the  above  definitions,  we  can  now  define  an  approximate  problem  as  follows: 
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Problem  (Pft):  Find  ^ e such  that 

(3.9)  J(v"h>  = inf  J(i|'.)  ' 


° *h  = tv,h  £ Sh!  - M'  ^h(0)  = fJ(0)'  30(3  *■  (1)  = f2  (0)  ). 


b)  J(t|/  ) = max  |g*  (nk)  - p(0)[f?(nk)  - f*  (nk)  ] 
n n_W0#  * * * ' WT  1 1 

- pto  ^^(0)  - p(°)  ^znh  (0)1. 

We  then  take  as  our  approxiamtion  to  u (t)  at  t = nk  the  function 

^,-r* 

(3-U)  U = W*  + Z + (1-x)  f * (nk)  + xf  * (nk)  . 

n n n 1 2 

We  now  show  how  Problem  (P^)  can  be  solved  by  linear  programming. 

Let  h be  such  that  H = 1/h  e TL  and 
H-l 

(3.12)  ij/h  = l C.$i  + f*  (0)  $Q  + f*(0)*H  , 


0 . x <_  (i-l)h, 

x/h  - (i-1)  , (i-l)h  < x < ih  , 

(3.13)  4>i  = / 

] 1+i  - x/h  , ih  <_  x <_  (i+l)h, 

0 , x > (i+l)h  . 

v*  — 

Note  that  the  constraint  1^1^  £ M is  equivalent  to  | | <_  M,  i = 1,...,H-1  and 

|fj(0)|  < M , i = 1,2.  Then,  by  linearity, 

H-  -r*  H-l  <0  f * (0)  $ +f  * (0)  $ . 

(3.14)  — Zn  (0)  = I C.  ±.  *>)  + ^ Z^  H(0)  - £ Znr  (0) 


i 3x  n ' ' 3x  n 


Hence  Problem  (P  ) can  be  written  as: 
A 


Find  <?  = (C^,...,C  ) minimizing  A subject  to  the  constraints 

a)  -M  < C.  < « , i = 1 , . . . ,H-1  , and 

(3.15)  b)  -A  < g*  (nk)  - p(0)[f*(nk)  - f* (nk) ] - p (0)  ■£-  W*  (0) 

2 1 9x  n 


H-l  4>.  f*(0)*„  +tt(0>*  , . 

-P(°)  (I  C iz  (0)  tj-z^1  H(0)  -^zr  (0)1 


n - No nt 


4.  Main  Results. 


In  order  to  derive  our  main  result,  we  shall  need  several  lemnas  about  the  regularity 
and  approximation  of  the  solutions  of  (2.1)- (2.2)  and  (2.3)  and  the  stability  of  problem 
(3.6).  We  shall  now  state  these  lemmas.  The  first  lenma  relates  the  smoothness  of  the 


solution  of  (2.3)  to  its  initial  snoothness  and  can  be  found  in  [2]. 


Lenma  2 : Let  be  the  solution  of  problem  (2.3).  Then  for  0 < tQ  < t and  s>0. 


there  exists  a constant  C,  such  that 

4 


(4.1) 


llzV.tJII  < C t"s/2|| 
s — 4 0 


Lemma  3 : (c.f.  Wheeler  [15]  ) . Let  w be  the  solution  of  (2.1)- (2.2)  and  {W  } be  its 


approximation  given  by  (3.6).  If  w e L (0,T;W 


,4 , 0 


3w  2 4 

e L (0,T;  H ) , and 


,3 
3 w 


— g-  € L2(0,T;H°),  then  there  are  constants  C.  and  k > 0 such  that  for  all  0 < k < k , 
a* J 4 0 0 


3t 


we  have  for  n = 0,...,NT, 


(4.2)  | w ( • ,nk)  - W I <c{h?[l|wll  . + II  II  _ 

n to  — 41  °°  4. 00  3 1 2 4 

1 L (0,T;W  ' ) IT  ( 0 , T ; H ) 


2„  3^w  „ 1 

+ 3 2 0 

3tJ  L (0,T;H  ) 


We  remark  that  sufficient  conditions  for  w to  have  the  regularity  required  by 
lenma  3 are  that 


a)  F e L2  (0 ,T ; H°) 


(4.3)  b) 


F(x,0)  = 0 for  x = 0 and  x = 1,  and 

qF(x.O)  + Ft(x,0)  = 0 , for  x = 0,1, 


. 3 , 3F(x,0) 

c)  3^  (p  ~ 


where  F is  given  in  (2.2).  These  are  natural  compatibility  assumptions  on  the  data  whose 
satisfactions  are  assumed  with  the  hypothesis  of  the  existence  of  a sufficiently  smooth 
classical  solution  of  the  problem  (1.1).  When  (4.3)  is  satisfied,  we  have  by  standard 


a priori  estimates  that  for  n = 0#...fN(j, 


(4.4) 


I w (•  ,nk) 


W I < C-{h?  + k2} 
n'«  — 5 1 


whe 


re  C,.  depends  only  upon  the  data  f^.f^,  p*  5 • Using  the  inverse  properties  (3.4) 


satisfied  by  the  subspace  S , we  easily  obtain  the  following  result  by  a standard  technique. 


Corollary  (3.1).  There  exists  a constant  C such  that  for  n = 0,...,N  , 


(4.5) 


| w ( • ,nk)  - W | < r^-  [lh  + k2] 

n 1 1 — h 1 


Our  next  lemma  allows  us  to  obtain  L -estimates  for  the  error  in  approximation  of 


(2.3)  by  (3.7)  for  times  bounded  away  from  t = 0 (i.e.  0 < t <_  t)  . 
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w 


Lr-nma  4:  u . f.  bdk  r , ft.  al  [1]).  Let  Z'  be  the  solution  of  (2.3)  and  (Z  j the 

— n 

approximat  i :•  (3.7).  Then  there  exists  constants  = C^(t^)  such  that  for 

n = N . . . . ,N  , 

0 T 

(4.6)  z^(*,nk)  - C7(h2  + K2  111  v»ll  Q . 

2 

Again , using  our  inverse  assumptions  on  S,  , we  obtain 

h 

Corollary  (4.1).  There  exists  a constant  such  that  for  n = 


| z v ( ■ ,nk)  - z’*l  < -A  [h2  + k2]lli>|l  . 

n 1 — h 0 


We  shall  also  require  a special  stability  result  for  the  approximation  scheme  (3.6). 
Since  the  proof  of  the  following  lemma  is  quite  technical,  we  shall  defer  it  until  after 

the  proof  of  our  main  result. 

Lemma  5 . Let  ) be  the  solution  of  (3.6).  If  for  some  constant  x Q > 0 , we  restrict 

k and  h^  such  that  kh^  _<  ^ , then  there  is  a constant  Cg  such  that 

<4-8)  iwJi  »_c9IIf"  „ o 

' L (0,T;H°). 

The  last  lemma  which  we  shall  state  gives  an  a priori  estimate  for  the  linear  pro- 
gramming problem  defined  in  (3.15).  Again  the  proof  of  this  lemma  will  be  deferred  until 

after  the  proof  of  our  main  result. 


Lemma  6.  Let  ^ 


be  the  solution  of  the  linear  programming  problem  (3.15).  Then 


there  exists  a constant  c such  that 


(4.9)  max  j< 

n=N  , . . . ,N 
O'  T 


g*(nk)-p(0)  [f*(nk)-f*(nk)]-p(0)-L-  W^(0)-p(0)^  Z^h  (0)  | = J(v>h) 


<_  C10U0  + h + [h^  + k2 ]/h^  + [h2  + k2]/hh 

We  are  now  in  a position  to  state  our  major  result  and  prove  it  using  Leninas  1-6. 
Theorem  1-  Let  u be  the  solution  of  Problem  (P)  and  {U^}  be  the  approximation  defined 
by  Problem  (Pft)  and  (3.11).  Suppose  that  hypotheses  (H1)-(H6)  are  satisfied,  that  F 
(defined  by  (2.2))  satisfies  the  regularity  conditons  (4.3)  and  that  if  = u(x,0)  satisfies 

noli ; : for  some  constant  > 0 . If  the  mesh  sizes  k,h,  and  are  chosen  to 

2 

satisfy  k = h = h^,  then  there  exists  a constant  C : which  is  independent  of  k such 

that  for  n = NQ>...,N  , 


|u(-,nk)  - U | < C [e  + kf 

n l , 00  — 12  0 


i 


where  eQ  and  y are  the  constants  defined  by  (H5)  and  (2.7)  respectively. 
Proof:  In  the  reformulation  of  problem  (P) , we  wrote  the  solution  as 

(4.11)  u = u^r  = w+z^r  + (l-x)f^(t)  + xf^tt)  . 

Fran  (3.11)  we  have 

-r* 

(4.12)  U = W*  + Z n + (1-x)  f * (nk)  + xf* (nk) . 

n n n 1 2 

Using  the  triangle  inequality,  we  obtain 


u(*,nk)  - U < w(*,nk)  - W + W - W* 

1 n 1 1 — 1 n i,“  n n'l,«° 

(4. 3.3)  f -r* 

+ |z*~r(-,nk)  - Znh  lx  „ + 2^0*)  - f * (nk)  | + 2 | f2  (nk ) - f*  (nk)  | . 


The  first  term  on  the  right  of  (4.13)  is  bounded  using  (4.5)  as  follows 
(4.14) 

Fran  (H5) , we  see  that 


| w ( • ,nk)  - W | < C [h*  + k2]/h  . 

n 1 ,OT  — 6 1 1 


(4.15) 


| f ^ (nk)  -f  * (nk)  | + 1 (nk) -f  * (nk)  | <_  2eQ 


Since  k = h^  by  hypothesis,  we  can  use  Lerrma  5,  (2.2)  and  (H5)  to  bound  the  second  term 
on  the  right  side  of  (4.13).  We  obtain 


(4.16) 


W - W*  < C II F— F*ll 

n n 1 >“  - 9 L-(0,TfH°) 


^C£o 


In  order  to  treat  the  third  term  on  the  right  of  (4.13)  we  use  the  triangle  inequality, 


(4.17) 


|z*-rf,nk)  - 

w.-r* 

£ | zf  r ( • ,nk)  - z (-fhk)^  + [z 


( • ,nk)  - Z 


■f.-r* 
h i 


Using  (3.3.b),  (3. 10. a),  and  (4.7)  we  have  for  n = Ng,...,N  , 


^h_r* 


( • ,nk)  - Z 


y>  -r* 
n 


1,-i  T ^ + k2)llVr*ll0 


1 1 


(4.18) 


(h2  + k2  ] { I y?  I + |r*|  } 
— h h «°  1 '■» 

< C[h2  + k2)/h  . 


Combining  the  above  estimates,  we  obtain 


(4.19) 


I u ( • ,nk)  - U I,  < C{e  + [h?  + k2]/h,  + [h2  + k2]/h) 

n l,00  — 0 1 1 


+ | -i  r (*,nk)  - z h ( • ,nk) 


1 ,°° 
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In  order  to  estimate  the  last  term  on  the  right  of  (4.19),  we  shall  use  the  continuous 


dependence  result  from  (2.7).  We  obtain 


(4.20)  \z  r ( • ,nk)  - z h (*,nk)|  < C_||s5*~r  (0,  •)  - * h (0,*)llY 

1,®  — 1 x x ' [t  ^ ,T]  * 

for  any  function  y(t)  e W ' [t^,T],  we  now  define  J^y  to  be  the  piecewise  linear  inter- 
polate of  y on  the  time  mesh  of  width  k . Then  we  obtain  (using  properties  of  a linear 


interpolate ) 


llzx'r(0'-)-ZIh  [t0,T] 


t4'21>  lllZx  r(0''>  ' ^lcZx~r  <0'  *,n  tt0,T]  + "Jkzx  h (0'°  - zxh  (0'-)ll[t  ,T] 


if  — IT  * 

+ n=W?“.,M  l<'r(0,nk)  - z h (0,nk)  | . 

0 T 


By  standard  results  in  approximation  theory,  the  first  two  terms  on  the  right  of  (4.2) 


bounded  by 


Ck2["4"t(0'-)"[t0,T,  +llzir  «>'•>«  [VT]> 


Next,  differentiating  equation  (2. 3. a)  and  using  Lemma  2 and  the  Sobolev  lerrma,  we  obtain 


for  any  ip  e H , 


llzxtt(0'->llro.Ti  icl^(-,t)li 


" xtt'“'  [0,T]  - 

(4.23) 

1 Cll  i|i  llQ  , for  0 < tQ  < t < T . 

Combining  the  above  estimates  and  using  (3.3)  and  (3.10.a)  to  see  that  ll^-r)l  Q < 2M  and 
II ,'’h“r*llo  i 2M  ' we  9et 

2 — r * 

1 C k + max  |z*_r(0,nk)  - z h (0,nk)  | 

(4.24)  n=N0,...,NT  x X 

- c k + max  |^-  zf  r(0,nk)  - j-  Z h (0)  | 

"=N0 Nt  3x  3x  n ' 


3 „Vr*. 


v d —r  * 

3 h , 


n=No NT 


4 


Using  (3.3),  (3.10),  and  (4.7)  we  estimate  the  last  term  on  the  right  of  (4.24)  by 

u-r"  c0 


n=N0 nt 


i;  zlh  r (0)  ' “ (0'nk)i  ^ f [h‘ + k‘>»*h-r*ii0 


(4.25) 


x 

. .2. 


£ Clh  + k ]/h  . 


Next,  using  (2.4)  and  the  triangle  inequality,  we  see  that  the  second  term  on  the  right  of 
(4.24)  is  bounded  by 


n=N  . . . . , N 
0 T 


3 z*'r(0,nk)  - ^ Znh  (0)  | 


3x 


n=N. , . . . , N_ 
0 T 


l2^-  It  (0'nk)  - tVnk)-fl<nk>)  - tZnh  (0)l 


— n=N N {FT5T  |g(nk)  - gMnk)|  + 'If  (0'nk)  - ^TWn(0)l 


max 

”0"  ’ * '*’T 


+ ' Jx  Wn(0)  ’37Vn(0)l  + lf2(nk)'f2(nk)  I + I f 1 ("*>  ~f  { <nk)  | } 


n=N  , . . . , N 
0 T 


|g*(nk)  - p(0)  (0)  - p(0)  [f  • (nk) -ff  (nk)  ] 


-p<°>  £ 


= term^  + term^  . 

Using  (H5),  (4.5),  and  (4.16)  we  can  bound  the  first  term  on  the  right  of  (4.26).  We 
obtain 
(4.27) 


term1  + [h^  + k2]/!^}. 


We  note  frcm  (3.9)  and  (3.10)  that  the  second  term  on  the  right  of  (4.26)  is  just  J(^) 
from  our  Problem  (Pft)  . We  then  use  Lemma  6 to  obtain  the  a priori  estimate 


(4.28) 


term2  = J(v>h)/p(0) 

-cio{£o  + h + thi  + k2]/hi  + th2  + *2]/h  }- 


Then  combining  (4.19),  (4.20),  and  (4.24)- (4.28)  , we  see  that  for  n = NQ,...,NT 

(4.29)  |u(*  ,nk)  - U |.  < C{  t + h + (h?  + k2]/h,  + [h2  + k2]/h}Y  . 

n l ,°°  — u 1 1 
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4 


Then  choosing  h k and  we  obtain 

(4.  V))  ,u(* ,nk)  - U 1 , < Ct  e -f  k}Y 

n 1 ,*  — 0 

tor  n N which  was  to  be  proved.  H 

Next,  assume  tnat  the  linear  programming  problem  described  in  Section  3 is  solved  to 

within  the  tolerance 

J(V 

for  some  ^ 0.  Replacing  the  estimate  (4.28)  by  the  above  inequality,  we  obtain  the 

following  error  estimate. 


Corollary  Tl:  Assume  the  hypotheses  of  Theorem  1 are  satisfied.  If  = V c*  $ 

h . i i 
1=1 

the  solution  of  t he  linear  programming  problem  defined  by  ( 3 . 15 )f  satisfies  J(<£*)  < , 

then  for  n we  have  for  some  constant  > 0 , 

'u(*'nk)  - UJl,-  -C13U0  + k + °1,Y- 

We  shall  finally  give  proofs  for  Lemmas  5 and  6 which  were  stated  previously. 
Proof  of  Lemma  b : To  prove  this  lemma  we  will  need  to  make  use  of  results  from  elliptic 
regularity  theory,  spectral  theory  in  Hilbert  spaces,  and  the  theory  of  interpolation 
spaces.  We  shall  assume  the  reader  is  familiar  with  these  concepts,  since  to  provide 
detailed  expl anat ions  would  unduly  lengthen  the  proof.  In  order  to  simplify  the  exposi- 


tion, we  first  introduce  some  additional  notation.  Let  Q be  the  solution  operator  for 
the  two  point  boundary  value  problem 

d)  '77(pTx)+qy=f'  0 < x < 1 , 

(4.  U)  b)  y(0)  * 0 , 
c)  ytl)  = 0 . 

( if).  Let  Q be  the  solution  operator  for  the  Galerkin  approximation 

hl 

of  (4.  31).  Then  y = Q f is  defined  by 
hi  hi 


aty.  ,vw  > 
hi  hi 


(f,v  ) , V e S . 
hl  hl  hl 


Now,  set  L « Q (i.e.,  the  inverse  of  Q on  S ).  For  normed  spaces  X and  Y, 
hl  hl  hl  hl 

let  l (X,Y)  denote  the  space  of  linear  operators  from  X to  Y and  II  *11 . 


X (X  ,Y)  the 


(operator)  norm  in  this  space.  We  shall  first  establish  that 

(4.33)  !w  I.  < C Ill7/8  W II  . . 

1 n 1 1,*  — n 0 


4 


We  note  that 


|w  |.  = | a l w | 

1 n 1,®  n l,a 


(4.34) 


i I \-V\  Vl.-  + l2\  WJl.» 

I I A.  K w - K 2 K - + K Q K w. 


hi  hi  n 


hi  hi  n i<“ 


w - 5 L w 

^ n n !,<■ 


+ Q L W I 

' n'l,«° 

1 4 

where  is  the  interpolation  operator  mapping  HQ  -*• 

4 1 * 

S.  given  in  (3.4),  we  see  that 

hl 

I Q.  L W - I Q L W I 
hl  hl  n hl  1 n 1»“ 


By  the  inverse  properties  of 


(4.35) 


-3/2 

1 c h IIQ.  \ w - I Q L w || 

1 n n 0 

u-3/2, 


< C h,  Mia  L W -01  W II  + II Q L,  W - I Q L W II  }. 

— 1 h^  h^  n hj  n 0 h^  n h M n 0 

Using  standard  properties  of  the  Galerkin  approximation,  we  note  that  the  first  term  on  the 

right  of  (4.35)  can  be  estimated  as  follows: 

(4.36)  h'3/2  lire.  - Q)L.  W II  < C h~3/2h?/4ll2  L W II,,,  . 

1 h^  * h^  n 0 — 1 l hin  7/4 

Also,  using  (3. 5. a),  we  obtain 


h"3/2|lQ  L.  W - X.  Q L W II  < C h"3/2h2/4  II Q L W II  _ . , . 
1 “hj^n  1^  h n 0 — 1 1 * h^  n 7/4 


(4.37) 

Then,  using  (3.5.b)  and  the  Sobolev  lemma,  we  see  that 

1 1.  Q L.  W - Q L W | < |X.  Q L W I,  + |q  L W |. 

1 h^  h^  n hi  n'l,”  — h h^  n'l,"  h^  n 1,“ 

(4.  38 ) < (C3  ♦ 1)  ICL^  W Jlf. 

-C  "2  \ Wn"  7/4  • 

We  then  collect  terms  in  the  above  inequalities  and  use  a standard  a priori  estimate  for 

-1/4  7/4 

elliptic  problems  (i.e.  Q e.  £ (H  , H ' ) (10)  ) to  obtain 


| w | < Cll  Q L W II  ... 

'n'l,®—  * h.  n 7/4 


(4.39) 


^Cll\  Wn 11 -1/4  ' 

where  the  last  norm  is  defined  below  in  (4.41).  Then,  in  order  to  obtain  (4.33)  we  need 
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only  show  that 


(4.40) 


II L W II  < C||Ly8  W II n 

hx  n -1/4  - hx  n 0 


2 4 2 4 

Let  P:L  S.  be  the  L projection  operator  into  S 

hl  1 


Note  that  L,  is  self- 
hi 


,1/8 


adjoint,  so  is  also  self-adjoint.  Then  we  see  that 

1 (L  W / v) 

h n 

II L W ||  . = sup  ij-ji 

1 v £ C 1/4 


SUP„ 

v £ C„ 


(L.  W ,Pv) 
h^  n 


II  vll 


(4.41) 


1/4 


= sup 

oo 

V £ C„ 


(Lhf  V Pv) 


inj 


1/4 


sup- 

1 v e C„ 


r i[6  pio 


II  vll 


1/4 


II  L?:/2  Pvll  = [a(Pv,Pv)]1/2  < dl  Pvll 
0 — 1 


Next,  since  we  have  that 

(4.42) 
and 

(4.43) 

we  can  use  an  interpolation  theorem  due  to  Heinz  (9)  (compare  also  with  [10)  ) to  establish 
that 

(4.44) 


iy  Pvll0  lllPvllo  • 


li  l^78  pvii„  < dl  pvll 


o - 


1/4 


Then , since  approximation  properties  of  S yield 

hl 


11  ^ 1/4  1 dlv,lx/4  - 


(4.45) 

we  can  combine  (4.41),  (4.44)  and  (4.45)  to  obtain 


(4.46)  II L VMI  < dl£-  Wnll0 


7/8 

1 " *'  ’ \ 

Now,  to  establish  (4.8)  we  solve  a set  of  first  order  difference  equations  arising  from 


(3.6)  to  obtain 


n-1 

W = k l 
n . *•- 

3=0 


(4.47) 


[1  + I*1  p F 


(•'  [n  - 7 - J1KM  . 


Then  we  see  tliat 


<[*  PFC.In  - i-  JR,II0 

n-1 


(4.48) 


♦k  I II  \/e  <11  + jU  1 -1 1 1 - k ] } j [ I + 7 L ] _1II 

1=1  hl  2 "l  2 "l  2 hl 


£ (H°,H°> 


• II  PF  ( • , [n  - j - j ] kll  Q . 


For  the  first  term  on  the  right  of  (4.48),  we  obtain 

"k  \[e  11  + i\'~l  PF(*'In  - r,k)llo 

<27/8k1/8i,,|,7/V/9  [i  + |t  r1 

~ ^ n.  o u 


(4.49) 


2 -h  ' " n n " PF  < tn  - r-]k)ll 

2 hl  £{H°,H°)  2 0 


7/8] 


< 27/8k1/8  sup  ( TT7  )llPF(*,  (n  - i-]k)|l 


0<X<°° 


<_  2 II  PF  ( - , In  - 2-1  k ) II  Q (for  k <_2). 


• —2 

Since  it  is  known  that  the  eigenvalues  A,  of  L.  satisfy  0 < A < 8 h , the  mesh 

i h 1—01 

-2  1 

ratio  restriction  kh  < t contained  in  the  statement  of  Lenina  5 implies  that 


(4.50) 


Ai  I 6oVk  • 


From  [1],  we  know  that  there  exist  constants  C and  K such  that  for  j = 1,2,...  , 

Vo 

and  0 < X _<  — - — , 


(4.50) 


nil 
1 + 1 x 


iv 


-r  i xj 
0 2 J 


“Cx  7 /8 

Then,  since  clearly  ex  is  bounded  for  all  x > 0 , we  use  (4.50)  to  see  that 


7/8 


-1 


V {tI  + 2\1  [I-  2 V}JtI  + 2 V1 


-1„ 


(4.51) 


hl  2 hl  2 hl  £ (H°,H°) 


sup 

0<X< 


boto 


1-ix 


1 ♦ ih 


>7/8 


1 + I x 


■C0  2‘\7/8 

0<  X<— 

— k 


sup 

0<X< 


_ k . , k 7/8 

"C0  2X]  ( 2A]) 


®oTo 


(jk/2) 


7/8 


1 C12(jk/2) 


-7/8  . 
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r 


Thus  combining  tfrie  above  estimates,  we  have 


7 /8  1 

II,  W II  < 211  PF  ( - , In  - -]k)ll  + llPFll 
i '.  ^ n 0 — 2 0 


0.  k J 


n-1  C 


12 


L (0,T;H  ) j = l (jk/2) 


>7/8 


< II  PFll 


12  fT  1 


L (0 


, T ; H° ) 12  + 7/8  I 


(4.52) 


2?/8  0 t?/8 


8C, 


dt] 


< II  PFll 


t2  + t1/8 
L~(0,T:H°>  27/8 


<_  c IlFll 

l“(0,T:H°) 

Finally,  combining  (4.39),  (4.46)  and  (4.52)  proves  Leirnna  5 . H 

Proof  of  lemma  6:  From  (H5),  (3.9),  (3.10),  (4.5)  and  (4.16)  we  see  that  for  all  ip  c K 

h h 

J(^h>  1 J(<!h) 

]p 

= nax  |g*  (nk)-p(O)  [f*  (nk)-f*  (nk)l-p(O)-^- W*  (0)-p(0)r~  Z h (0)  I 

n=N N 2 1 3x  n ^ 3x  n 1 

0'  ' T 

1 max  { | g*  (nk)  -g  (nk)  | + p (0)  | f * (nk)  -f  (nk)  | + p (0)  | f * (nk) -f,  (nk)  I 
n=N0 \ 

+ P(0,l  ^Wn(0)  ' 37Wn(0)l  + p(0)  Wn(°)  ' £ 


+ P(0)|^rZnh  (0)  * 37  z^"r(0,nk)|} 


Now 


(4.54) 


42  3 *h_r*  3 

< C[eq  + (hx  + k J/l^]  + max  p (0)  1 3^"  zn  (0>  ~ z^~r  (0,nk) 

n = N0 \ 


, 0, -r*  . 

\~-  Z h (0)  - ~~  z^”r  (0,nk) 
1 3x  n 3x 


, \l>.-r*  . 'lv_r*  3 'I'. -V>tr-r* 

1 |—  Z_  (0)  - — z (0,nk)  | + |gr  z (0,nk)  | . 


— 1 3x  n 


3x 


3x 


From  (4.7),  (3.3),  and  (3.10)  we  see  that  for  n = Ng,...,N  , 


(4.55) 


t|/  -r*  ip.-r* 

id  _ h 3 h 

— Z (0)  - — z 

3x  n 3x 


< C[h2  + k2l/h. 


(0 ,nk) | < [h2  + k2]IU-r*l 

n h 0 
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Next  using  the  Sobolev  lenma,  Lemma  2,  and  (3.3.c),  we  obtain  for  n - N(|,...,N  , and  for 
any  \ \ - 

* -v+r-r*  \k  -v>+r-r* 

| r — z (0,nk)  | <_  II  z II 


(4.56) 

< Cll  ik  -i^tr-r*ll 
— n u 

i c[t0  t IUh  - *tQ\  . 

Choose 

tj/,  such 
h 

that 

a) 

i k (ih)  = tf  (ih)  , i = 1, . . 

h 

. . , H-l 

(4.57) 

b) 

H-h(0)  = fj(0)  , 

c) 

H>h(l)  = f*(0)  . 

Then  il  e K.  and  i k is  almost  equal  to  I.'f,  the  piecewise  linear  interpolate  of  •f  . 
h n h h 

Then  from  (3.5), 

(4.58)  <.  11  \J»h-Xhv»H  o + IIJ'hv’^’ll0 

< II  ik -I.  ^11  . + Ch  Ml,  . 

— h h 0 1 

But  we  see  that 

[ 0 ' h 1 x 1 (H-l)h  , 

(4.59)  = < (1  - ^)  (f  J (0) -f  x (0)  ] , 0 <_  x <_  h , 

[ (H-l)  ] [f*  (0)  -f  (0)  ] , (H-l)h  < x < 1 , 

^ h 2 2 — — 

so  we  have 

(4.60)  11  *h  - Vo  iC  E0  • 

Combining  (4.  53)  - (4 . 60)  , we  obtain 

(4.61)  J(¥>  ) < C(e„  + h + (ft  + k2 ) /h  + (h2  + k2)/h]  . * 

h — o i l 
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